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ABSTRACT 

If Jupiter and the Sun both formed directly from the same well-mixed proto- 
solar nebula, then their atmospheric compositions should be similar. However, 
direct sampling of Jupiter's troposphere indicates that it is e nriched in element s 



20081 ). 



such as C, N, S, Ar, Kr, and Xe by 2-6 x relative to the Sun (IWong et al\ 
Most existing models to explain this enrichment require an extremely cold proto- 
solar nebula which allows these heavy elements to condense, and cannot easily 
explain the observed variations between these species. We find that Jupiter's 
atmospheric composition may be explained if the Solar System's disk heteroge- 
neously accretes small amounts of enriched material such as supernova ejecta 
from the interstellar medium during Jupiter's formation. Our results are similar 
to, but substantially larger than, isotopic anomalies in terrestrial material that 
indicate the Solar System formed from multiple distinct reservoirs of materia l 



Trinquier et al. 



2003). 



simultaneously with one or more nearby supernovas {e.g 
Such temporal and spatial heterogeneities could have been common at the time 
of the Solar System's formation, rather than the cloud having a purely well-mixed 
'solar nebula' composition. 



1. Introduction 

The Solar System's composition reflects that of the initial cloud core from which it 
formed. On small bodies such as the terrestrial planets, subsequent thermal and chemical 
processes have altered the composition. However, for the Solar System's two largest bodies 
- the Sun and Jupiter (mass 1 Mj 0.001 Mq) - the atmospheric composition is relatively 
stable against change since the time of formation. In particular, the relative atmospheric 
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abundances of gases carbon (C), nitrogen (N), sulfur (S), argon (Ar), krypton (Kr), and 
xenon (Xe) are believed to be fixed at the time of formation 4.5 billion years ago, and 
not affected substantially by subsequent evolution. The noble gases are stable against loss 
due to Jupiter's high escape velocity, and stable against sinking by condensation because 
they don't combine with other atomic species. Elements C, N, and S combine into larger 
molecules such as CH4, NH3, and H2S, but are believed to remain well mixed in Jupiter's 
troposphere without settling to the interior. 



Jupiter's composition at 5-20 bars was 



directly measured b y the Galileo descent probe's 



2003 



hereafter AMN03) showed 



mass spectrometer in 1995. Initial results ( lAtreya et al\ 
that Jupiter was enh a nced by ~ 3 ± Ix relative to the then-current solar abundances of 



Anders and Grevessd (119891 . hereafter AG89). More recent work has better determined the 
Solar abundances, leading to revised Jupiter: Solar abundances. Somewhat surprisingly, 
these new result s actually increase d the scatter between the elements, to the broader 



range of 4 ± 2 x (IWong et al. 



20081 . hereafter WLA08). The largest ch ange occurred in the 



abundance of Jovian Ar, which more than doubled to 5.4 ± 1.1 x solar fiGrevesse et al. 



20071 . 



hereafter GAS 07) . A subsequent measurement of solar Ar reduces it again to 2.45 ± 0.7 x 



fcodders 



20081 ) (hereafter Lod08). All enrichments in this paper are in terms of number 



density relative to hydrogen, normalized to the Sun; i.e., {ni/nu) / {ni/nH)s 



2. Previous models 



Several previous models have been proposed to address Jupiter's 'metallicity problem. 
All are based on condensation of the volatile elements into ic es, which can then b e 



concentrated rela tive to hydrogen. The 'amorphous ice model' (jOwen and Encrenaz 



Owen et al. 



2006 



I999I ) proposes that Jupiter's atmosphere was formed from a combination 



of solar nebula material and 'solar composition icy planetesimals' (SCIPs). The SCIPs 
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were made of amorphous water ice in which the other volatile species condensed; these 
SCIPs had the same composition as the Sun, except were depleted in H. Jupiter then 
formed from a linear combination of these two components: when SCIPs were introduced 
into Jupiter from the protoplanetary disk, they were heated and released their volatiles 
into Jupiter's atmosphere. (This two-component model is distinct from t he formation o 



Jupiter's co re, which is a separate issue.) The 'clathrate hydrate model' (iHersant et al. 



2008 



20041 ) uses a similar approach, but proposes that the volatiles are trapped within 



clathrate hydrates, rather than amorphous water ice. In this model, species condense out as 
the nebula slowly cools to their clathration temperature. The process can work at slightly 
warmer temperatures than the amorphous ice model. The condensation order of hydrates 
leads naturally to the Jovian abundance. Both models have the appeal of simplicity and 
directness. 

These models share a common problem, that the nebula must be very cold for volatiles 
to condense. The condensation temperatures required are 20-35 K for th e noble gases 



onto amorphous ice or into solids, and 35-55 K for the clathrate hydrates (ilro et al. 



2003 



Hersant et al. 



2004j ). Th e clathrate models a ssume very cold nebula temperatures of <20 K 



at 5 AU at several Myr ( iHersant et al. 



20041 ). These temp eratures are qu i te low compared 



with both observations and models of disks. For instance, 



20-50 K at 50 AU, and much warmer inward, a nd 



100 AU. Models by 



work by 



Dullemond and Dominik 



the model of ( IJonkheid et al. 



Watson et al. 



Dutrey et al. 



( 120071 ) measures 



(|2003[ ) measures 50 K at 



Chiang and GoldreichI ( Il997l ) predict 75 K at 5 AU and 55 K at 10 AU; 



20051 ) assumes 90 K at 5 AU and 30 K at 50 AU; and 



20071 ) predicts 50 K at 100 AU. Large uncertainties exist in 
both observational and theoretical studies of disks, and it remains to be established that 
temperatures can indeed reach the low temperatures needed by these models. The low disk 
temperatures required are also dangerously close to the 20-50 K gas and dust temperatures 
of the molecular cloud itself, which set an absolute lower limit for the disk temperatures 
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( iJohnstone and Ballyl l2006l : iJohnstone et all 120061 : iBally et al\ Il99ll ) . The molecular cloud 
temperatures are highest in the largest clouds and those with massive stars, which are the 



regions most likely to be where the Solar System formed (iHester et al 



20041 ) ■ 



Several models address the temperature problem by forming solids not at their current 



location, but further out in the nebula {e.g., 10-50 AU) where temp eratures were 



Radial migration subsequently brought them to 5 AU. For instance, 



Alibert et al. 



ower. 



fl2005h 



explores the possibility that Jupite r formed in its enti r ety at 10-15 AU via the clathrate 



model, and then migrated inward. iGuillot and Huesol (120061) use photo-evaporation and 



viscous migration of the entire nebula to explain Jupiter's enrichment. In this work, UV 
flux from internal and/or external stars heats and preferentially removes hydrogen from 
the inner proto-planetary disk. These planetesimals then migrate inward where they heat 
and go into forming Jupiter. This model is appealing, and consistent with the idea that 
the young solar system niay have experienced the effects of nearby massive stars {e.g., 



Tachibana and Huss 



2003 



Throop et al. 



200 If ). 



Bringing material from much further out, where the nebula remains cold, is appealing 
but there is a strong dynamical constraint on h ow much material from the c oldest regions 



of the disk could be incorporated into the disk. iLevison and Duncan 



fll997[ ) showed that 



less than 0.5% of comets leaving the Kuiper Belt eventually impact Jupiter. To acquire 
the amount of contamination seen here would require that the Kuiper Belt have a mass 
on the order 2500 Earth masses. Models of Neptune's outward radial migr ation to 30 A U 
place firm upper limits on the Kuiper belt mass of just 20 Earth masses (iGomes et al. 



2004 



Hahn and Malhotra 



19991 ). so Jupiter likely received no more than 0.1 Earth masses 



of material from the Kuiper belt or beyond. 

Even if the proper temperatures are somehow reached, additional problems exist for 
the clathrate model. Laboratory results at very low pressures are sparse and have not 
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shown that clathrat ion rates at the very low nebula pressures are sufficiently fast to operate 



on Myr timescales f lHersant et al 



20081 ). And, in order for ice to trap other volatiles as 



clathrates, its surface must remain exposed ('microscopic' grains or cracks) for several Myr 



until the nebula cools ( jHersant et al. 



20041 ) ■ This is inconsistent with models showing that 



grains grow rapidly; pl anetesimals up to 100 km are thought to be able to form within 



1 Myr even at 30 AU (iWeidenschilling 



19971 ). If for any of these reasons clathration is 



inefficient, then Jupiter's core mass could grow much higher than observed due to high 0:H 
ratios required. Finally, this existence of such small grains would work to insulate the disk 
and keep it warm, rather than allow it to cool as required. Thus, the solar system's bodies 
may simply have formed too large, too quickly, to clathrate sufficient volatiles from the disk 
atmosphere. 

All the models discussed do a reasonable job of fitting Jupiter's composition. The 
SCIP models predict a uniform enhancement of 3x in all species, fitting the original Galileo 
data very well. However, no SCIPs have ever been identified in the Solar System: comets, 
for instance, ar e depleted substantially in N and are of decidedly non-SCIP composition 



( liro et al. 



20031 ). The clathrate models result in too much S (by a factor 2x). All of the 



models attempt to fit the roughly 3 ± Ix enhancement of the origina l Galileo values, and 



do not try to fit the newer 4 ± 2 x values; 



Owen and EncrenazI ( l2006l ) dismisses the revised 



values as being possibly due to systematic errors. None of the models fit the high Ar value 
of 5.4 ± 1.1 X of GAS07, though the 2.5 ± 0.5 value of Lod08 clearly fits more easily. Ar has 
in fact the coldest condensation temperature (~35 K into clathrates, and 20 K as a solid) 
of all the species measured, so its high GAS07 value would be a particular challenge to all 
the condensation models. 



Finally, 



Lodderd (120041 ) proposed an entirely different model: that Jupiter formed from 



carbon-rich (rather than ice-rich) planetesimals. However, this scenario depends on the 
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assumption that the Gahleo probe sampled a typical region of Jupiter's atmosphere, and 
most work supports the opposite view that the probe hit an anomalously dry spot. 



3. Our model: A contaminated molecular cloud 

We propose an entirely different solution to the problem. Rather than forming Jupiter 
and the Sun from identical 'Solar nebula' material and invoking condensation or transport 
within the nebula to modify Jupiter's composition, we propose that the composition of 
Jupiter and the Sun differ because of intrinsic temporal and spatial variations in the 
Solar nebula composition as the interstellar medium (ISM) is polluted by massive stellar 
winds and supernovae (SNs). We show that this model can explain Jupiter's chemical 
composition, easily fits into environmental formation scenarios, and is consistent with other 
heterogeneities in both our Solar System and distant star clusters. 

In the scenario that we propose, multiple stages of star formation occurred within 
a giant molecular cloud (Fig. [1]). The cloud was of average size, 10-20 pc, and within 
the GMC lay several pc-scale molecular clouds. Stars, including the Sun with its disk, 
began to form in clusters within these clouds. The Sun's orbit through the cluster took 
it on a long path several pc across. Before any nearby 0/B stars turned on, the cluster 
environment remained cool and dark for several Myr. During this time the Sun passed 
through t he ISM and gravitationa lly swept up material onto its disk by Bondi-Hoyle 



accretion (IThroop and Bally 



20081 ). 



If the Sun formed alone, then the composition of the ISM might be uniform. But in our 
proposed scenario, other clusters nearby formed a few Myr earlier, some with higher-mass 
stars, and the ISM surrounding these clusters soon began to be polluted in heavy elements 
produced by these massive stars. Nucleosynthesis in these stars enriched their ejecta by 
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~100x or more relative to Solar H. The Sun's evolving orbit took it through these more 
polluted regions of the ISM. Material was accreted onto both the Sun and its disk, but the 
disk's large cross-section and low mass made it easier to pollute than the Sun. The disk's 
metallicity slowly increased, and Jupiter's core and atmosphere formed from the disk. The 
disk dispersed within 5-10 Myr, around the same time as the local ISM dispersed and the 
cluster spread apart, ejecting the Sun as an unbound field star. Jupiter's final composition 
thus represented largely the same material as the Sun, but with the late accretion of 
polluted material reflecting in Jupiter's enriched composition today. 

Our proposed model reflects much of the curren t understanding of star formation in 



large clusters and molecular clouds {e.g. iBallyl 120081 ). Stars do not form in isolation, the 



ISM is not of uniform composition, stars travel on long orbits through their young clusters, 
and material from the ISM can be accreted onto young stars and disks in the several 
Myr after they form. All of these are newly appreciated processes that have not been 
incorporated into existing models of Solar System formation. We flnd Jupiter's metallicity 
to be one natural consequence of these processes, and we explore here the issues involved 
with it. In the following sections we examine the detailed chemical constraints on Jupiter's 
composition from such a model (S]), and the timing and spatial parameters required (^. 
In §5.31 we describe several other cases where stellar pollution occurs. 



4. Chemical constraints 

In this section we investigate the detailed chemical enrichment from this 'polluted 
accretion' scenario and whether it can explain Jupiter's current composition. 

We start by assuming that Jupiter's atmosphere is composed of a mixture of three 
distinct components. The majority of the mass is made of Solar material, whose composition 
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is well determined by the measurements of GAS07 and AG89. The remainder is the 
small amounts of pollution that come from massive stellar winds and/or supernovae. The 
compositions of these ejecta are distinct 'fingerprints' of the stars, determined mostly 
by their initial mass and composition. In general, the winds of these stars are enriched 
primarily in light elements (C, N, O, Ne), while the SN ejecta several Myr later contain 
heavier species (S, Ar, Kr, Xe). As one example, the demise of a star with an initial mass 
of 20 Mq can produce 5 of O, 0.5 M,^ of C. 0.2 of Si, 0.001 M© of ^SAl, and 10"^ 



Mq of Fe in its SN phase ( iWoosley and Heger 



20071 . hereafter WH07). In our models, we 



use the grid of winds and ejecta computed by WH07. These models span the mass range 
12-40 Mq, and consider a variety of nucleosynthesis rate coefficients and explosion energy 
parameters, for a total of 66 models. 

We denote the elemental abundances in the Solar, wind, and SN components as n©, 
?7,„, and risn, the Jupiter abundance is rij. It is then possible to calculate model Jupiter 
compositions by using 

'^J.i = /© nQ,i + /w ^w,i + /sn '^sn.i; (1) 

where i is the species and /© + /sn + /w = 1- We use the present-day solar abundance, 
but the difference between this and the Sun's primordial composition is insignificant for our 
purposes. 

Using Eq. [H we found combinations of solar material and ejecta that would produce 
Jupiter's measured composition. We searched all 66^ = 4356 possible combinations of the 
66 wind and 66 SN models of WH07, coupled with the single solar abundance. For each 
trial, we computed coefficients that best fit Jupiter. Our routine attempted to fit only the 
well measured stable species (C, N, S, Ar, Kr, Xe), and computed results for both these 
and the remaining elements (He, O, Ne, P). 

Our best fit ('Model A') is shown in Figures[2H3l This model finds Jupiter's composition 
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to be well described by 87% solar nebula, 9% stellar winds from a 40 Mq star (WH07's 
model s40a28A), and 4% supernova ejecta from a 20 Mq star (WHOT's s20a37n). The 
total contamination is 13% (^.e., 0.13 Mj). The fit is excellent at matching the observed 
quantities of C, S, Ar, and Kr, and the lower limit for O. The largest deviation is for N, 
where we are slightly below the error bar. The wind predominantly supplies C, N, and O 
while the SN supplies the remaining species. Both stars have high enough mass (and thus 
short enough lifetimes) that they can form and explode within the 10 Myr timeframe of 
GMCs. We assume the GAS07 and WLA08 values for the Solar and Jovian composition. 
This model is shown here because it is the best fit; many other combination of different 
mass SN and wind ejecta provided much worse fits. 

A variant of this fit ('Model A2') is shown in Figures HHS] In this case we have used 
the latest LodOS value for Ar, instead of that of GAS07. This model is fit with 78% solar 
nebula, 8% winds from a 40 Mq star (WH07's s40a28A), and 14% SN ejecta from a 15 Mq 
star (WH07's sl5a34c). The only difference to the fit is the Ar abundance. Because 
lower-mass SNs are less efficient at heavy-element nucleosynthesis, the fit here requires a 
substantially larger SN contribution than does Model A (14% vs. 4%). The 15 Mq star has 
a lifetime of ~11 Myr, putting it on the upper end of individual cloud lifetimes but within 
the timescale of large regions like Orion. 

Finally, a third fit is shown in Figures [6H7] ('Model B'). This model differs in that we 
have used the abundances of AG89 and AMN03 for rigoi and njup. Although the newer 
abundances are probably preferred, using the old ones gives a test of the robustness of our 
fits. Also, the nucleosynthetic yields of WH07 start with the AG89 solar abundances, so in 
a sense this fit is more self-consistent, even though it is based on slightly older data. Model 
B requires about 6% total contamination, less than half that required by Model A. Model 
B is comprised of 94% solar, 4% stellar winds from a 40 Mq star (s40a28A), and 1.5% SN 
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ejecta from a 25 Mq star (s25a41d). 

All three of our models fit the data well. In all cases the vast majority of Jupiter's 
mass (78%-95%) comes from the Solar nebula, in agreement with standard models. And 
in each, a combination of winds and SN - which one would expect in a realistic cluster - 
works better than any single component by itself. All three fits are slightly low in N and 
Xe; the remaining species are fit very well. 

The individual 'fingerprints' of the SNs can be seen in Figures [71 |5l and [31 In all the 
models, the lighter species (C, N, O) are produced in the winds, while the heavier elements 
come from the SNs and are highly dependent on the SN mass. A large difference between 
the three models is the SN C:N ratio, which is several times higher in Model B than the 
others (compare red curves in [SI and [3]) . The great deal of carbon ejected by the Model B 
SN allows the total contamination in this model to be about half that in Model A. 

Detailed yields from our three models are listed in Table [H This table lists additional 
species measured by Galileo, but which are not in equilibrium at the entry site and 
thus not expected to fit: He, Ne, and P. Helium and Ne are believed to combine into 
He-Ne 'raindrops' which sink t o Jupiter's interior, and cannot be used as a constraint 



( jRoulston and Stevenson 



19951 ). Helium itself is produced by the Sun so its primordial 
Solar abundance cannot be directly measured. O was measured at Jupiter but its value 
is believed to be anomalously low due to the probe's dry entry site. Encouragingly, our 
models predict primordial Jovian abundances for O, Ne, and P similar to those of the other 
species. The predicted values for O are in the range 2.6-2.9, similar to global O values 
determined spectroscopically (WLA08). 

Several changes could improve the quality of our fits. First, we have used a simple 
fitting method, assuming contamination by only one wind and one supernova. In realistic 
star-forming regions such as Orion there are several dozen stars within a few pc all above 
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8 Mq that will explode as SN; using multiple stars will increase the ease of fitting Jupiter's 
composition. Second, the SN ejecta models we use are quantized in relatively large mass 
bins (5 Mq), and all use certain common assumptions for stellar and expl osion parameters 



he S N ejecta yields are very model dependent; for instance, the models of 



Young and Fryer 



( 120071 ) vary in abundance for individual species by 50% or more from those of WH07 for 



stars of similar mas s. The WHO? yields ignore stellar rotation, which may be important 



(IHirschi et al. 



20051 ). Thus, a better understanding of SN yields could help (or hurt) our fits. 
Finally, we have not included physical processes in the ISM or disk such as sedimentation, 
concentration, or condensation, even though these could be important. For instance, of the 
noble gases our poorest fit is for Xe, which is low by 50% and at the edge of the 1 a error 
bar. However, Xe condenses at a substantially warmer temperature than any of the other 
noble gases (at ~55 K into clathrates, it is the easiest to condense), so if our model were to 
include condensation explicitly, it would operate in the direction to correct this deficiency. 
For now, however, we have intentionally chosen to keep our models simple to demonstrate 
that good fits are possible even with a very limited set of parameters. 



5. Stellar Pollution into the ISM 

Now that we have shown that Jupiter can be matched chemically, we investigate the 
requirements on the environment to support such contamination. 

As massive stars in a cluster evolve, they perform nucleosynthesis and create heavy 
elements that pollute the ISM. These are given off in stellar winds (during the stellar 
lifetime) or SNs (at the end-of-life for stars with M > 8 Mq). After being injected into 
the ISM, these highly enriched ejecta are incorporated into the next generation of stars. 
The highest-mass stars are the shortest lived: for instance, the 40 Mq stars used in our 
fits explode in less than 5 Myr, allowing for a full generation of star formation within the 
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10 Myr timescales of planet formation and embedded clusters. 



5.1. Contamination by Stellar Winds 

Stars spend the majority of their lives on the main sequence. During this stage the 
stellar winds are weak enough that they do not pollute the cluster; for instance, the current 
solar wind is Mq ^ 2 x 10~^^Mq yr~^ with v 400 km s"^. Red supergiant (RSG) 
stars, however, have slow massive winds that can easily pollute the cluster. RSGs are 
normal post-main-sequence stars of mass 5-15 Mq that have cooled dramatically after their 
main-sequence phase (T > 10,000 K to T < 3000 K), and have heavy and slow winds of 

10~^-10~^ Mq y r~^ and v as low as 10 km s~^ - ten order of magnitude more loss than 



the present-day Sun (IKnapp and Woodhams 



1993h. The RSG ph ase lasts for 5-10% of the 



stellar lifetime, or typically a few 10^-10^ yr (iSchaller et a/.iil992iV For stars > 25 Mg^, ove r 



19961). 



half the original stellar mass can be lost during the RSG phase (iGarcia-Segura et al 
The wind velocity of 10 km s~^ is only slightly higher than the stellar and gas velocities, so 
this material can readily mix with the local ISM. The abundances n„ that we consider in 
§ m are predominantly from the RSG phase, because stars of less than 20 Mq evolve too 
slowly (> 10 Myr) to enter the RSG phase within typical cluster lifetimes. Our Model A 
and B fits use winds from stars of 40 Mq with lifetimes ~ 5 Myr, which easily fit the timing 
constraints. These high-mass stars are seen in clusters of > a few thousand; an example 
is 6^ Ori C in the Orion Trapezium core. Smaller 8 Mq stars are produced in clusters of 
N > a few hundred. 
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5.2. Contamination by Supernova Ejecta 



After passing through the RSG phase, stars with masses > 8 Mq usually end their 
lives as SNs. The SN explosion gives off 1-5 Mq of metal-enriched materi al at speeds of 



2,000- 10,000 km s ^. Even at such high speeds, numerical simulations by 



Quellette et al. 



( 120071 ) showed that disks at 1 pc are very resilient to nearby explosions. However, although 
the disk will absorb most of the intercepted solids such as ^°Fe, most of the SN's gas is 
deflected by the disk. For a 20 Mq explosion at 1 pc, the disk will absorb just 10"'* Mj of 
enriched gas, insufficient for the process described here. 

But, there are at least two mechanisms by which ejecta may cool and then be accreted 
onto the disk. First, as the ejecta spreads, it mixes with the ISM until it slows and cools, 
and this cooler ejecta can be accreted more easily. In order to slow to 10 km s"^, 1 Mq 
of SN ejecta must mix with roughly 1000 Mq of ISM {i.e., a 1000:1 mixing ratio). The 
SN ejecta can be quite clumpy, with high mixing ratios in some regions and low in others. 
Rather than mixing uniformly, these clumps appear to be slowed as a unit and preserve 
their density, much like a baseball is slowed in the air without fully mixing. Observations 



of pc-scale ejecta from various SNs s how highly c 



maintained for a year after explosion ( iFesen et al. 



umpy knots on AU scale s or greater 



2007 



Wang and Hu 



19941 ). Observations 



of emission lines even 700 years after the explosion continue to find evidence for density 



contrasts o 



Fesen et al. 



100- 200 in dense metal-rich bubbles within the ejecta (jWilliams et al. 



2008 



Myr timescales (IFreyer et al. 



2006h. Additiona 



wor k shows that some h eterogeneity persists on even longer 



2006 



Kroeger et al. 



20061 ). Thus, diluting to a 1000:1 mixing 



ratio is needed to slow the material, there will likely be local regions where the mixing ratio 
could be 100:1 or lower, roughly the amount of SN pollution required by our 'Model B' 
case. The CMC itself is generally robust against disruption by SNs, as evidenced b y the 



Orion region, which has seen dozens of supernovae in the past 12 Myr (IBally 



2008|), yet 



- 15 - 



remains largely a cool GMC 



The second mechanism involves AU-scale 'droplets' of ejecta that can travel from 
a supernova exploding as a 'champagne fountain,' through the ISM, and land cool 
and intact in molecular clouds many kpc distant. This mechanism is believed to be 
responsible for observed compositional heter ogeneities ('abundance discrepancy factors ' 



or ADFs) of 50%-lOx in molecu. 



2007 



Tsamis and Pequignot 



2005 



ar clouds (iMesa-Del 



Tenorio-Tagk 



lado et al. 



2008 



Stasinska et al 



19961 ) ■ These ADFs are seen in cold. 



dense, hydrogen-depleted gaseous regions in the middle of larger molecular clouds, on 
scales of 10-100 AU. They are commonly seen in star-forming clouds including Orion 
and 30 Doradus. The droplets have mixing lifetimes of several Myr, consistent with the 
timescales needed for our Jupiter model. In this way, dense concentrations of gaseous 
metals produced by the SN are transported to and preserve d in the nebula, where they are 



then passed on to stars and planets that form in the region (jStasinska et al. 



20071 ). 



5.3. Accretion of the ISM onto disks 



Once enriched material has mixed wi th the ISM by eit 



accretion onto the disk is straightforward. 



ler of these two methods, 



Throop and Bally 



( 120081 ) found that the average 



ISM-to-disk accretion rate for disks in young clusters was 10 ^ Mq yr ^, or ~1 MMSN per 
Myr. The total accretion needed in our Model B case is 0.1 Mj = 10"^ MMSN = 10"^ M©, 
or the amount delivered in ~ 10,000 years. Since accretion may be maintained for 5 Myr, 
this provides more than sufficient mass delivery. Due to the disk's large cross-section 
and its low mass, accretion affects the d isk but makes only negligible effect on the stellar 



composition (IMoeckel and Throop 



20091 ). 
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6. Heterogeneity in the Solar System and star- forming regions 



Our model reflects the growing body of evidence that the Solar System did not form 
from a homogeneous cloud in an isolated environment, but rather from a heterogeneous 
nebula where interactions with its environment played a major role in shaping its evolution. 
We briefly discuss here four examples of large-scale heterogeneity: two in the terrestrial 
planets, one in star-forming regions, and one on galactic scales. 

The terrestrial bodies have been modifled by scores of chemical and physical processes 
since their formation. Most of these processes act equally on all isotopes of the same 
element, so isotopic differences are not expected in samples formed from a well- mixed 
nebula. However, isotopic variations have indeed been measured in samples from the 
Earth, Mars, and numerous asteroids. Isotopic differences have been measured for species 



including Ba, Cr, S, Ti, Z i, Mb, O, and more (ITrinquier et al. 



2006 



2007 



Ranen and Jacobsen 



Dauphas et al. 



2002 



and references therein). These anomalies are small (up to 1% for 



^^Cr/^^Cr, and as small as ppm for some others) but are in disputable. UV photoch emistry 



has been invoked to explain the origin of the O variations (ILyons and Young 



20051), but 



the remaining heterogeneities have defled explanation by known fractionation processes. 
Instead, they are consistently thought to be of nucleosynthetic origin, resulting from the 
i ncomplete mixing of eject a from multiple SNs in the material of the young solar nebula 



(iRanen and Jacobsen 



2OO2I). 



2008 



Trinquier et al. 



2007 : 



Ranen and Jacobsen 



2006 



Dauphas et al. 



In addition to the stable isotopes, variations have been found in short-lived radioactive 
isotopes, al l owing some constraint on the timing of the Solar System's early heterogeneity. 



Krot et al. 



(120081 ) examined the heterogeneity of ^^Al isotopes within Ca-Al inclusions 
(CAls) in carbonaceous chondrites. They found differences of > 100 x in primordial ^^Al 
abundances, leading to the conclusion that there are at least two populations of CAls: 
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some which were formed in the presence of ^^Al, and some which were not. This requires 
either spatial or temporal variations in the Solar System birth environment. Their preferred 
interpretation is that the ^^Al-free CAIs were formed early (possibly during initial collapse 
of the Solar System), followed by late injection of ^^Al from a nearby massive star. They 
suggested that this massive star could be the same one that injected ^°Fe several Myr later 
after an SN explosion. We caveat this point with mention that a recent study of CAIs 



withi n a single ordinary chondrite did not reproduce the ^^Al variations ( jVilleneuve et al. 



20091 ): it is possible there are regions that were heterogeneous and regions that were not. 



Orion, the nearest massive star-forming region, shows evidence for 'polluted accretion' 
amongst individual young stars. Observations of the metallicity of 29 B, F, and G stars 
within Orion have f ound abundance variations up to 4 x between stars of the same age in 



the same subgroup (ICunha et al. 



2000 



1998 



Cunha and Lambert 



19941 ). These variations 



are spatially correlated, suggesting that recent supernovae have contaminated distinct 
regions of the cluster. The abundance variations are seen in O and Si (which are produced 
by massive stars and type II SNs), but not in Fe, C, and N (which are produced to a far 
lesser degree, and thus would not be expected to contaminate the stars). The authors 
proposed that ejecta from recent SNs have been accreted onto these stars, causing their 
observed metallicity. 

Relative to the galactic value, the entire Solar System exhibits a system-wide excess of 
30% in ^^O/^^O, and the solar value is higher than nearly every molecular cloud or YSO 
within 10 kpc. This enhancement has been proposed to be due to accretion of supernova 
eiecta and massive stella r winds immediately prior to the proto-Solar nebula's formation 



( lYoung et al. 



2009 



20081 ). 



These four examples describe heterogeneity broadly similar to what we propose for 
Jupiter. The contamination at Jupiter is many times higher than the isotopic anomalies in 
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the terrestrial bodies. However, the terrestrial anomalies were formed earlier in the solar 
system's history, over a shorter period, and in a warmer and better-mixed region. Jupiter 
formed later and in the outer solar system, where chemical pollution can more easily 
be acquired from the ISM, and mixing times are far longer. While the small terrestrial 
anomalies were likely to have been inherited from incomplete mixing of the initial solar 
nebula, Jupiter's enhancement could be caused by post- collapse contamination of the disk, 
from Bondi-Hoyle accretion onto the disk over several Myr and several pc, allowing for far 
more introduced pollution across the entire molecular cloud. The amount of heterogeneity 
seen at Jupiter is comparable or smaller than that seen in both star-forming regions (ADFs) 
and young stars themselves. Because the physical processes of star formation are universal, 
if the Sun formed in a dense cluster, then the solar system and Jupiter could have naturally 
inherited these heterogeneities. 



7. Discussion 

The model we propose here is a departure from existing models for the solution to 
Jupiter's metallicity problem. Although historically most Solar System formation models 
have assumed a homogeneous 'Solar nebula' composition, our model explicitly assumes the 
opposite. We incorporate the fact that the composition of the solar nebula can change 
spatially and temporally. These changes occur during the Solar System's first 5-10 Myr, 
as the Sun is traveling through its birth cloud, experiencing the environmental effects of 
other stars. The inclusion of this contamination reflects the latest understandings of the 
environmental processes affecting star formation. 

Our model has three general advantages over the existing amorphous ice and clathrate 
models that we describe in §[2j First, it relaxes the very strict low-temperature requirements 
for the formation of Jupiter's solids. Second, it allows for Jupiter to have formed at its 
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present location without migration. Third, it explains the fact that different elements are 
enriched by different amounts {i.e., not a uniform 3x). 

Our model is not without problems. Our fits for N are at the edge (or a little beyond) 
the la level. Nitrogen is poorly measured in the Galileo probe data, but is a major species 
and an important constraint. We could increase the contribution from stellar winds (which 
supply much of the N), but this would increase the amount of C beyond that observed. 
Low-mass stars of 4-8 Mq produce large quantities of N sufficient to solve the problem, 
but these stars have lifetimes of 30 Myr or more before they enter the RSG phase. The 
other specie our models are all low in, Xe, has a particularly high condensation temperature 
(making it easy to condense), and this may in part explain its abundance. 

Second, the amount of SN contamination required is on the high end of what can be 
supplied by nearby SNs using direct injection of ejecta into a nearby GMC Additional 
work is necessary to understand the mixing and concentration mechanisms of this ejecta. 
However, if the 'droplet' model of SN ejecta is correct, then concentrated ejecta can be 
slowed and cooled over kpc distances, which would solve this problem. 

Most work suggests that Saturn, Uranus, and Neptune accreted their atmospheres 
from the disk in much the same way Jupiter did. The metallicity enrichments of the outer 
planets exceed Jupiter's - Uranus and Neptune have 30 times the Solar C:H ratio, for 
instance - suggestin g that they were formed at least in part by the traditional chemical 



condensation series ( iLewis 



19721 1. The colder nebula and slower evolution make volatile 
condensation easier at greater distances. However, 'accreted pollution' of the disk from the 
ISM would affect these planets as well. Because of their different positions in the disk and 
their different formation times, our model cannot predict the enrichment they may receive 
from pollution, except that it could be similar to Jupiter's. (The outer disk, with its larger 
cross-section, might be more easily contaminated, but this requires more study.) More 
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detailed analysis of this awaits better knowledge of their bulk atmospheric compositions, 
which are largely unknown today. Jupiter's oxygen composition will be probed by Juno 
in 2017 and will provide a discriminant between many models including our own (which 
predicts a global 0:H ratio of ~ 2.6 x solar) and that of Lodders 2004 (which assumes 
~ 0.35 X solar). 

Jupiter's atmosphere may have formed through a variety of processes, of which polluted 
accretion could be only one. For instance, N is among the most difficult species to condense 
in the various ice models, as it requires a temperature < 30 K. Nitrogen, however, is easily 
supplied in the cool, slow stellar winds that we study here. It is possible that N was 
supplied by these winds, and other species were delivered in part by icy planetesimals from 
a relatively warm (> 50 K) disk. 

An advantage of the SN model is that it provides a mechanism for matching not only 
the chemical abundances which we study here, but also isotopic differences. It has already 
been shown that SNs may explain the isotopic differences seen in the Solar System's rocky 
bodies (§IS]); additional work on the isotopic differences between the Sun and Jupiter would 
be a natural extension of the present work. This path is particularly rip e for future study 



given the recent So lar composition results from the Genesis mission (e.q.. iWiens et al. 



2007 



Mabry et al. 



2003). 



The general process of accretion from ISM disk — > planet that we describe here has 
broader applications for the formation of extrasolar planets. A strong observed correlation 
exists between th e metallicity of the host star and the existence of extrasolar planets {e.g., 
Udry et a/.l 120071 ). We predict that the correlation between stellar and planet metallicities 
will depend on the formation environment. Stars that form in dense regions, in the presence 
of massive stars, will have the highest disk metallicities and thus the greatest chance of 
forming planets. 
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Fig. 1. — The proposed 'polluted accretion' scenario for Jupiter's atmosphere. In this model, 
the Sun and its disk form in a low-metallicity molecular cloud (left). The Sun's orbit takes 
it through other regions (right) of higher metallicity, polluted by massive stellar winds and 
supernovae. Ongoing Bondi-Hoyle accretion from the ISM delivers this enriched material to 
the disk, where it is incorporated into Jupiter's atmosphere. 
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Fig. 2. — Elemental abundances of Jupiter from Galileo probe (grey), and fits for our 'Model 
A' case (red). The y axis plots the elemental number abundances relative to the Solar 
abundances, normalized to hydrogen. The model consists of a linear combination of 87% 
solar composition, 9% from stellar winds from a 40 Mq star, and 4% ejecta from an SN 
of original mass 20 Mq. O is a lower limit because the Galileo probe entered Jupiter at 
a cloud-free location believed to be anomalously dry. Solar and Jupiter compositions are 
revised 2007 values (GAS07, WLA08). 
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Fig. 3. — Sources of individual species in our 'Model A' case. The model is a linear combi- 
nation of elemental abundances from the Sun (green curve), stellar winds (blue curve), and 
an SN (red curve). Each line is normalized to H at 1.0. The stellar winds produce much of 
the C and N, while the supernova supplies most of the remaining species. 
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Fig. 4. — Elemental abundances of Jupiter from Galileo probe (grey), and fits for our 'Model 
A2' case (red). The model consists of a linear combination of 78% solar composition, 8% 
from stellar winds from a 40 Mq star, and 14% ejecta from an SN of original mass 20 Mq. 
Abundances are the same as Model A, except the Lod08 argon value is used. 
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Fig. 5. — Sources of individual species in our 'Model A2' case. Same as Figure [3l but using 
Lod08 argon value. 
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Fig. 6. — Our 'Model B' fits. Same as Figured but assuming 1989 solar composition data 
(AMN03, AG89). The coefficients are Solar (94%), stellar winds from 40 Mq star (4%), and 
a supernova from a star of original mass 25 Mq (1.5%). The fractional contamination in 
this model is 6%, less than half that in Model A. The low predicted value for Xe might be 
explained by its particularly high condensation temperature (see text). The broad similarity 
of the results to the 'Model A' case shows that our model is robust against small changes to 
knowledge of the composition of the Sun and Jupiter. 
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Fig. 7. — Individual components of our 'Model B' fits. Same as Figure [3], but using 1989 
solar abundances (AMN03, AG89). 
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4% 




WLAOS 


8% 


14% 




WLAOS+LodOS 


4% 


1% 




AMN03 


H / 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


1.0 


He 


1.0 


2.1 


2.3 


1.2 


0.8+0.0 


2.1 


1.9 


1.2 


0.8+0.0 


2.1 


3.7 


1.1 


O.S +0.0 


C V 


1.0 


33.9 


9.5 


4.3 


4.3 ±1 


33.9 


6.4 


4.3 


4.3 ± 1 


33.9 


36.6 


2.9 


2.9 ±0.5 


N / 


1.0 


6.7 


6.5 


1.7 


4.9 ±1.9 


6.7 


4.9 


2.0 


4.9 ± 1.9 


6.7 


10.3 


1.4 


3.2 ±1.4 


O 


1.0 


6.3 


31.3 


2.6 


0.3 ±0.2 


6.3 


11.3 


2.9 


0.3 ±0.2 


6.3 


9S.I 


2.6 


0.3+0.2 


Nc 


1 il 
















1 1 _ _ 1 1 






J 1 


i 1 _ 1 1 ( 1 


P 


1.0 


1.6 


170.0 


7.5 


0.8 ±0.0 


1.6 


17.6 


3.4 


0.8 ±0.0 


1.6 


139.9 


3.1 


0.8 ±0 


S / 


1.0 


1.5 


71.0 


3.7 


2.9 ±0.7 


1.5 


13.4 


2.8 


2.9 ±0.7 


1.5 


106.9 


2.6 


2.5+0.15 


Ar / 


1.0 


1.5 


103.4 


4.9 


3.4 ±1.1 


1.5 


11.5 


2.5 


2.5 ±0.7 


1.5 


93.3 


2.4 


2.5+0.5 


Kr / 


1.0 


1.5 


16.6 


1.6 


2.2 ±0.4 


1.5 


92 


2.2 


2.2 ±0.4 


1.5 


118.4 


2.7 


2.7 ±0.5 


Xe V 


1.0 


1.5 


1.7 


1.1 


2.1 ±0.4 


1.5 


1.4 


1.1 


2.1 ±0.4 


1.5 


2.9 


1.1 


2.6 ±0.5 



Table 1: Details of model results, based on the listed ejecta yields from winds and SN 
compared with Jupiter. Table lists species observed by Galileo; abundances are in number 
density relative to H, normalized to Solar. The checked species are those which are stable at 
Jupiter and our model attempts to fit; we predict the unchecked species but do not attempt 
to fit them. The three models are based on different measurements for the Jupiter:Solar 
abundances. 
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